Continuum spin excitations in S=l one-dimensional antiferromagnet. 
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We present detailed measurement of the non-hydrodynamic part of the spin excitation spectrum 
in the model quasi-lD S=l antiferromagnet CsNiCl3 by inelastic magnetic neutron scattering. In 
the better part of the Brillouin zone the effect of the inter-chain coupling is negligible, and spin 
dynamics reflects that of a single Haldane chain. We find that at q < 0.6-7T this quantum spin system 
ceases to support coherent propagating excitation, which gradually turns into a continuum band of 
states, whose width grows with decreasing q. This finding is consistently verified under different 
resolution conditions, obtained with two different high-luminosity experimental setups. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.50.Ee 
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I. INTRODUCTION 

In striking neglect of the naive "common wisdom", 
Haldane ground state (GS) of the one-dimensional (ID) 
S=l Heisenberg antiferromagnet (HAFM) does 
not connect the Neel-ordered "spin solid" GS of the semi- 
classical 1 HAFM with the almost ordered "marginal 
liquid" state of the S=l/2 chain O]. Instead, it is a 
"quantum liquid" with finite correlation length and a 
gap in the spin excitation spectrum. Spectral weight 
of the spin fluctuations is concentrated in a long-lived 
massive triplet mode in the neighborhood of the Bril- 
louin zone (BZ) boundary q = n. Any remainder of 
the spectacular continuum observed at q w tt in the 
S=1/21D HAFM H is predicted to be extremely faint 
P,p[p^]. On the other hand, both nonlinear er-model, 
which is believed to be a valid description of the S=l 
HAFM chain in the long- wavelength limit, and the vari- 
ational treatment based on the Jordan- Wigner fermion- 
ization, predict that two-magnon continuum states are 
the lowest-energy excitations at q w and dominate 
the spectrum in the vicinity of the Brilloun zone cen- 
ter |l4],[L5]]. In the absence of magnon interaction the 
two-magnon continuum starts above a threshold energy 
£2m(q) = min{2 ■ e(n + q/2), e(ir + q) + Ah}, smallest 
at q = 0, where it is 2A#, twice the Haldane gap. 

Experimental observation of the continuum part of 
spin excitation spectrum in Haldane chain is a very chal- 
lenging problem. Primarily this is due to the rapid de- 
crease of the static structure factor S(q), which gives the 
energy-integrated intensity of the scattering cross-section 
by spin fluctuations, at small q ||[4J,0|. In single- mode 
approximation (SMA) S(q) ~ (1 — cosq) / E(q) vanishes 
~ q 2 as q — > 0. In fact, in the state of the art neutron 
scattering experiment JtJ authors found no sign of a con- 
tinuum in the model Haldane chain compound NENP 
down to q = 0.3-7T. This result has lead to a widespread 
belief that there is indeed no continuum at q > 0.37T, 
despite indications to the contrary by MC studies . 
However, upon careful examination the findings of the 
Rcf. |?J are not at all that prohibitive. To merely ob- 
serve magnetic scattering at q = 0.3ir authors employed 



very coarse resolution, with full width at half maximum 
(FWHM) of the instrument wavevector acceptance cov- 
ering about a quarter of the ID BZ. Their statement that 
measured magnetic intensity is consistent with the SMA 
cross-section, treated within the context of this broad 
resolution, only imposes an upper limit, and not a very 
stringent one, on the width of the continuum. Another 
obstacle to observing the continuum magnetic scatter- 
ing is a single-ion anisotropy in the spin Hamiltonian. 
As in the majority of Ni-organic chain compounds, it is 
quite large in NENP, and results in a significant splitting 
of spin fluctuations with different polarizations. Smeared 
by the resolution such splitting is hard to distinguish from 
continuum. 

Quasi- ID antiferromagnet CSN1CI3 is one of the most 
isotropic and best studied Haldane model compounds. 
It provided some of the earliest, although initially con- 
troversial, experimental evidence in favor of the Haldane 
conjecture B. In CsNiCi3 super-critical inter-chain ex- 
change coupling J' w 0.03J results in a long-range or- 
der below Tn w 4.8 K, but as temperature rises above 
Tjv, a gap opens in the spin excitation spectrum, and 
it quickly recovers properties of an isolated S—l HAFM 
chain Unfortunately, characterization of T = 

continuum is impossible at such elevated temperatures 
because of the substantial magnon thermal broadening. 
However, it was recently pointed out on the basis of 
the MF-RPA (mean field random phase approximation) 
analysis that inter-chain coupling modifies only the low- 
energy part of the excitation spectrum, and therefore 
even at T < TV dynamic spin response of CSMCI3 in 
the better part of the BZ around the top of the ID dis- 
persion is identical to that of an individual chain Jl6[ . 
Here we present neutron scattering measurements which 
illustrate this point and reveal the existence and the ex- 
tent of the excitation continuum. 



II. EXPERIMENT 

We studied large m rs 6.4 g sample of CSMCI3, com- 
posed of two single crystals, co-aligned to yield effective 
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mosaic < 1° . Sample was mounted on an Al plate in the 
standard "ILL orange" 70 mm cryostat with (h, h, I) zone 
in the scattering plane. CSMCI3 has hexagonal structure 
P63 / mmc with two equivalent ions per c spacing, so that 
Q = (h,k,l) in reciprocal lattice units (rlu) corresponds 
to gii = nl in the ID BZ of a chain. Sample had longer di- 
mension parallel to the hexagonal c-axis, and was in the 
transmission "Laue" geometry shown in Fig. 1. Mag- 
netic scattering was measured at the "base" temperature 
T = 1.5(2) K, the non-magnetic background (BG) was 
collected in identical scans at T = 150 K. 

Experiments were performed on SPINS 3-axis cold 
neutron spectrometer at NIST Center for Neutron Re- 
search, using two complementary setups illustrated in 
Figure 1 (a), (b). Monochromatic incident neutron beam 
was obtained by (002) reflection from vertically focused 
pyrolytic graphite (PG) monochromator, viewing the 
58 Ni neutron guide. While in setup (a) the "natural" 
beam collimation around the sample was « 100' — 540', 
it was restricted to 80' — 80' by Soller collimators in setup 
(b). Scattered neutron wavevector was analyzed by (002) 
reflection from an array of 11 (a) or 9 (b) independently 
rotating w 2 cm wide PG crystals through a 80' radial 
collimator (RC) onto rj 24.3 cm wide position sensitive 
detector (PSD) with 256 pixels. Calibration of the scat- 
tered neutron energy accepted by PSD pixel and its sen- 
sitivity was done by measuring clastic incoherent scatter- 
ing from standard sample at different incident energies. 




FIG. 1. Two high count rate setups with PSD used for 
measurement of the ID spin excitations in CsNiCl3 . (a) large 
flat PG analyzer gives about 10-fold increase in acceptance to 
the direction of the scattered neutrons, (b) 9 PG segments 
aligned to create a polychromatic image of the sample at the 
RC center similarly increase acceptance to the length of k/. 

Use of the high count rate setups which employ SPINS 
large area segmented PG analyzer and matching large 
PSD was key to the success of our measurements. Essen- 
tial element in such setup is a radial collimator in front of 
the PSD. It dramatically restricts acceptance of the PSD 
pixels to incoherent scattering from the analyzer and its 
support structures, which otherwise render phohibitively 
large background. RC requires that reflection from the 



analyzer produce a geometrical image of the sample at 
its curvature center. For each incident wavevector the 
direction and/or length of the scattered wavevector k/ 
are functions of the pixel position across PSD, as defined 
by geometry of the Bragg reflection at the analyzer, Fig. 
1. This results in a coupled scan with both energy trans- 
fer E and wavevector q in the sample reciprocal space 
varying across the PSD. However, q component along 
dispersive direction of interest (chain in our case) can be 
kept qu w const in such scan, if this direction is aligned 
perpendicular to the direction of the scan by appropri- 
ately choosing q±, as shown in Fig. 1. 

Setup (a). PSD central energy was fixed at — 
4.2 meV, and range Ef G [3.6,4.8] meV was covered 
in this measurement. Evidently, large flat analyzer sim- 
ply creates a (polychromatic) mirror image of the sam- 
ple. Radius of the SPINS RC is equal to the sample- 
detector flight path, optimizing this setup. For all PSD 
pixels (all sample scattering angles 29 s ) k/ component 
along the analyzer Bragg wavevector t a is constant, 
= Ta/2. Therefore, in the paraxial approximation con- 
stant is q\\ « < - k ' r '^- 4 - ) — 2^, while q± varies across PSD as 
a function of k f . 

Setup (b). PSD central energy was fixed at E^ = 
4.57 meV. Angle between the consecutive PG segments 
5(9a) ~ 35' -7- 40', and the corresponding difference in 
the reflected neutron energy, were chosen to create a 
polychromatic image of the sample at the RC center. 9 
analyzer segments provided full PSD coverage, reflect- 
ing energies in the range Ef G [4.03,5.13] meV. Chain 
direction was aligned perpendicular to kj, and there- 
fore, for given kj and 29 s , only q± varied across PSD, 
while qu = ki sin 29 s — const. Important feature of this 
setup is a possibility to restrict analyzer angular accep- 
tance by inserting a collimator after the sample. This al- 
lows measurement at smaller scattering angles, and also 
shapes/tightens the instrument resolution function, as 
shown by the resolution FWHM ellipses in Fig. 2. 

III. RESULTS AND DISCUSSION 

In Figure 2 we show contour plots of the measured 
spectral density of the magnetic scattering intensity, 
I(q,E)/ J I(q,E)dE, with the linear (q, independent 
background subtracted. It is evident in both panels that 
spectrum acquires finite width in energy at I < 0.5. 

Fig 2 is constructed from the raw data (integration is 
done via point by point summation), and is slightly dis- 
torted by the instrument resolution. In the single-mode 
part of the spectrum, at I > 0.6, measured line-shape is 
completely defined by the interplay of dispersion and res- 
olution. Although resolution volume is smaller in setup 
(b), the "focusing" effect (longer axis of the FWHM el- 
lipse is parallel to the dispersion in Fig. 2(a)) results in 
sharper peaks in setup (a). In principle, an opposite "de- 
focusing" effect is of concern for I < 0.5 measurements in 
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this setup, as it would result in quite significant broaden- 
ing even of a single-mode spectrum. However, careful ac- 
counting for the resolution shows that non-zero intrinsic 
width at I < 0.5 accounts for > 2/3 of the spectral width 
measured in setup (a). "De- focusing" is absent in setup 
(b), where the FWHM ellipse is approximately round, 
Fig 2(b). Spectacular agreement of the q\\ < 0.57T spec- 
trum measured in two setups shows explicitly, that ob- 
served cross-over from a single mode to continuum with 
significant intrinsic width is not a resolution effect. 
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FIG. 2. Contour plot of the spectral density of magnetic 
scattering reconstructed from the constant-^ scans, mea- 
sured in setups (a) and (b) of Fig. 1 correspondingly, via 
linear interpolation. Scale on the right shows variation with 
energy of the wavevector transfer perpendicular to the chain 
at I — 0.5, scale on the top - its variation with I at E — 3 meV. 
Ellipses are calculated half maximum contours of the instru- 
ment resolution function at E = 3 meV. Solid curve is the 
single-magnon dispersion discussed in text, dashed line shows 
the lowest energy of two non-interacting magnons with given 
total 9|| = irl, dotted line is e^ii) = 2. 49 J sin qu . 

Another important distinction between the data in Fig 
2(a) and (b), is that qu w const geometry imposes differ- 
ent conditions on q± in setups (a) and (b) correspond- 
ingly. Therefore, measurements are done along different 
trajectories in the sample 3D reciprocal space, as illus- 
trated by the right and top axes in Fig 2. The effect of 



the inter-chain coupling is manifested at low energies by 
difference in dispersion at I > 0.8. While in (a) h ~ 0.3 
is close to the 3D magnetic Bragg position at I » 1, 
for setup (b) it is close to the top of inter-chain disper- 
sion, h « 0. In agreement with MF-RPA analysis, ef- 
fect of the inter-chain coupling becomes insignificant at 
0.25 < I < 0.75. For comparison, curves in Fig 2 show 
magnon spectrum in a Haldane chain obtained in [15| , 
for J = 2.275 meV, Ah = 0.41 J , v = 2.49 J, a = v |§. 

Our results show explicitly and unambiguously, that 
single-mode excitation becomes unstable around the top 
of the dispersion band in a Haldane chain. Instead, a 
continuum excitation spectrum, whose width increases 
with decreasing q, is observed at q < 0.6-7T. In fact, very 
similar behavior was found in another quantum liquid 
- superfluid 4 He, where the "maxon" excitation turns 
into a broad continuum-like feature under pressure |p7j |, 
which suppresses the "roton" gap and drives the system 
towards crystallization, a quantum phase transition. 
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